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We present analytical expressions for the eigenstates and eigenvalues of electrons confined in 
a graphene monolayer in the presence of a disclination. The calculations are performed in the 
continuum limit approximation in the vicinity of the Dirac points, solving Dirac equation by freezing 
out the carrier radial motion. We include the effect of an external magnetic field and show the 
appearence of Aharonov-Bohm oscillation and find out the conditions of gapped and gapless states 
in the spectrum. We show that the gauge field due to a disclination lifts the orbital degeneracy 
originating from the existence of two valleys. The broken valley degeneracy has a clear signature on 
quantum oscillations and wave packet dynamics. 


I. INTRODUCTION 

The presence of topological defects in condensed mat¬ 
ter systems can modify thoroughly their physical proper¬ 
tied. From the point of view of the geometrical effects, 
topological defects can induce curvature and/or torsion 
in the geometrical background. Differential geometry can 
thus be used for the study of a variety of physical sys¬ 
tems, ranging from superfluid helium and its numerous 
spontaneous symmetry breakingd, to liquid crystald or 
mesoscopic physics 1 4 3. The prominent work of Katanaev 
and VolovictP has showed the efficiency of this geometric 
approach to investigate the propagation of elastic waves 
in the presence of topological defects such as screw dis¬ 
locations, edge dislocations or disclinations. 

In a recent paper, one of us showed how the presence of 
a screw dislocation can generate spin current in an elec¬ 
tron gad. A screw dislocation is a source of torsion. Tor¬ 
sion can also be built in directly in the physical structure, 
like it is the case in chiral molecules where chirality cou¬ 
ples to spin in the presence of spin-orbit interaction and 
can be used for spin selectivity applications, e.g. to po¬ 
larize an electron beanP or to exalt a spin orientation of 
a spin current propagating through the molecul^ 7 . Edge 
dislocations are also the source of torsion, and their role 
was studied in topological insulators 8 . 

It is desirable to consider in this context topological 
defects which generate curvature. Single disclinations 
are line sources of curvature and they correspond to the 
generation of a conical geometry via the Volterra con¬ 
struction. Our aim in the present paper is to investi¬ 
gate the properties of a graphene layer in which the crys¬ 
tal symmetry is locally modified (a hexagon is replaced 
by a pentagon or a square for example), giving rise to 
large scale properties analogous to the conical geome¬ 
try mentioned above. The desire to employ graphene 
for studying curvature effects is motivated by the sim¬ 
plicity of the Hamiltonian and the important potential 
applications of graphene in naoelectronics and possibly 
future quantum computing deviced. Carbon nanotubes 
are other famous elementary exemples of curvature ef¬ 
fects in a carbon sheets, but the case of disclinations is 


different due to the presence of a singularity which sits 
at the defect location. 

Recently, there have been many studies on quantum 
rings, in which the confinement of electrons with phase 
coherence of wave function gives ris e the Aharonov- 
Bohm 10 and Aharonov-Casher effect^ 112 . Quantum 
rings have been investigated both experimentally and 
theoretically in semiconductor deviced^ and also in 
graphene layers both in monolayers and bilayers 14 18 1 In 
this context it is important to study the curvature effects 
on electron properties e.g. charge currents or quantum 
oscillations in a graphene monolayer ring. 

Graphene also offers a way to probe quantum field 
theory in condensed matter systems, as its low energy 
spectrum is described by the Dirac Hamiltonian of mass¬ 
less fermions. The Dirac electrons in graphene occur in 
two degenerate families which correspond to the presence 
of two different valleys in the band structure - a phe¬ 
nomenon known as fermion doubling. Due to the valley 
degeneracy, in many cases, it is difficult to observe the in¬ 
trinsic physics of a single valley in experiments, because 
the two valleys have equal and opposite contributions to 
a measurable quantity. One way to break the valley de¬ 
generacy is the production of fictitious magnetic field in a 
single valley by a lattice defect. The field has an opposite 
sign in the other valley, so it lifts the degeneracy. The 
broken valley degeneracies leads to some interesting fea¬ 
tures on quantum oscillations and revival times of wave 
packets, which are the main concern of this paper. We 
show that the signature of the broken valley degeneracy 
is clearly visible in the persistent current and the wave 
packet characteristic oscillation times. 

Here, we consider a model of graphene with a wedge- 
dislocation, which can be understood from Volterra’s cut- 
and-glue constructions. Such conical defects are elemen¬ 
tary objects which are observed in graphene layers and 
in various carbon based structures 19 23 . Our calcula¬ 
tions are based on the continuum model where fictitious 
gauge fields due to the defects are coupled to the Dirac 
fermions 1 24 [ We solve this model by freezing out the 
carriers’radial motion, keeping into account only their 
angular motion in a effectively one-dimensional system 
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(namely a quantum ring). The full calculations includ¬ 
ing carriers radial motion will be presented elsewhere 25 . 

The paper is organized as follows: in Sec. [II] we present 
the theoretical model and numerical results for the energy 
spectrum, pseudo-spin polarization and charge current. 
In Sec. m the revivals and typical time scales of electrons 
wave function are presented. Sec. |IV| contains a summary 
of the main results and conclusions. 


II. THEORY 

In the low energy limit, the Hamiltonian of a graphene 
layer is given by 

H = v F [r z a x p x + (jyP y \ + A t z g z (1) 

where & = (cr x , cr y ,cr z ), r = (r x ,r y ,r z ) are Pauli matrices 
denoting the sublattice and valley degrees of freedom, 
respectively, and v F is the Fermi velocity. The finite 
mass term A, which might be caused by an interaction 
with the underlying sub strate or by spin-orbit coupling 
is also taken into account 26 2 ^. The Hamiltonian ([!]) acts 
on the four-component spinors 

= (i>A+(r) ipB+(r) ipA-(r) ip B -(r)) T - (2) 

The pseudospin indices A and B label the two sublattices 
of graphene and the valley indices + and — refer to the 
two inequivalent Dirac point K and K' in the Brillouin 
zone. In a two dimensional plane for a flat graphene, the 
angular boundary condition for a Dirac spinor as it goes 
a closed path is given by 

T(r, 0 = 2 t r) = e i7r(7zTz ^(r, 0 = 0) (3) 

The presence of defects like disclination modifies the 
angular boundary conditions in Eq.([3|. It is indeed 
known that deformation in the honeycomb lattice en¬ 
ter in the continuum limit description as fictitious gauge 
fields cou pled to the electron momentum as a vector 
potential 1 28 * 29 l We introduced a wedge disclination via 
the Volterra construction: first a sector of titt/3 is re¬ 
moved or added on the flat graphene lattice and this 
is then followed by the gluing the edges of the cut (see 
Fig. [l]). An index n defines the type of disclination in 
the honeycomb lattice. Keeping the symmetry of the 
lattice, a hexagon can now be replaced by pentagon 
(n = 1), square (n = 2), heptagon (n = —1) and oc¬ 
tagon (n = — 2) by the cut-and-glue procedure. 

From Eq.pk, the wave function is multiplied by the 
factor to keep the Hamiltonian covariance. So, 

when a sector is removed from the graphene plane, the 
wave function should transform according to 

(r, 9 = 2tt) = -e i2nl 1 ~( n / 6 )^/ 2 ^(r, 9 = 0 ) (4) 

When n is odd, an additional phase is required for the 
parallel transport of a Dirac fermion since when the 
spinor is transported around the apex, the phase mis¬ 
match due to the sublattice permutation (as for example 


a B site comes to a B site instead of A) must be com¬ 
pensated to preserve the complete wave function single¬ 
valued. Considering all the factors, the angu lar boundary 
condition for a general n now become^ 20 * 30 * 

= 27 r) = - e i2 ^~^ na y T y^^^ 1 ~^ n ^ 6 ^ azTz ^ 2 ^(r, 0 =(0) 

Note that for even n, the matrix r y does not play any 
role. Eq. © can be written as 

T(r, 0 = 2tt) = W n V(r, 0 = 0) (6) 

where 

W n = e iS r( azTz - 3a y Ty / 2 ) m ( 7 ) 

Where we rescale the angle 0 of the unfolded plane to the 
new angle 0 = 0/( 1 — f), where 6 is varying from 0 to 2tt. 
The general gauge transformation given in Eq. © for any 
type of disclination is defined by the index n = ±1, ±2 
defined above. For a local gauge transformation, Eq. © 
can be written as two singular gauge transformations 

T(r, 0 = 2tt) = U(</>)V n (0)*{r, 0 = 0), (8) 

where J7(0) = e L ^ (JzTz and V n (0) = e l ^ T y (J y. 



(9 = 0, 0 = 0 


= 2t r, 0 = 2 tt - ^ 


FIG. 1. Unfolded plane of lattice where a wedge of angle 
utt/3 is removed (n = 1 here). A potential (sketch on the left 
part of the ring) confines the electrons on a ring of radius p, 
avoiding the singularity at the origin. We rescaled the angle 
0 of the unfolded plane to the new angle 0 = 0/(1 — f). 


Let us now consider an additional external magnetic 
flux <F applied at the centre of the defect hole (a magnetic 
flux tube) corresponding to an azimuthal vector potential 


A = 


<!> 




(9) 


with £l n = (1 — ^) and <Fo = f, the magnetic flux quan¬ 
tum. We write the Hamiltonian in Eq. 0 in polar coor¬ 
dinate and apply the gauge transformation through U (0) 
and V n (0). The transformed Hamiltonian is denoted by 
H(r,9) = WVjHV n U, 


H(r,9) = hv F (k r - — )r 2 cr ;r 


+hv F (kg + 


$ 


+ 


f2„$ 0 r 4 tt n r 


-T z )(J y + At z ct110) 
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where k r = — and ko = Written in matrix form it is 


rVt n dO 


J 


H{r,0) = 


rA 

-i{d r + J) + ^r-(dt) + i§- + Ijt) 


~i(dr + £) - d e + i§~ 0 + * jO 
—rA 


( 11 ) 


where r = +1 and —1 correspond to the K and K' points 
respectively. The eigenstates of Eq. are two compo¬ 
nent spinors which, for each valley index, are given in 
polar coordinates by 

"> = '“’(,!$)) < 12 > 

where J is the total angular momentum and ^(r) and 
^(r) are linear combinations of Bessel functions due to 
the cylinder geometry. Now we will proceed the calcula¬ 
tions by freezing the radial motion. Considering a one¬ 
dimensional ring and simplifying the Schrodinger equa¬ 
tion by discarding the radial variation of electron wave- 
function is a correct approximation for a narrow ring for 
which the dynamics is frozen in the lowest radial mode. 
Excited radial modes have an energy much higher and 
can be neglected 31 ! when one studies low energy proper¬ 
ties. The radial part of the momentum is thus set to zero 
(P p = ^ 0 in matrix element^ 32 * 33 !). In 

the case of a ring of fixed radius r = p, this leads to the 
substitution -A — ^ in Eq. (10) and the eigenstate 


(12) becomes 


T r (r = p,0) = e 


jo ( $a(p) \ 

V*Mp)7 


(13) 


where &a(p) and ®b(p) are now fixed amplitudes. Due 
to cylindrical symmetry, the spinors (13) are eigenfunc¬ 
tions of the total angular momentum J z = L z + S z with 
S z = \h<J z and h(m+ is the eigenvalue of total angu¬ 
lar momentum (m G Z). People have often used different 


I 

notations for the wave functiorP^ 16 18 i like 


*t(p,0>® 


imO 


$A(p)e 

i<f> B (p)e i(m+1)e 


(14) 


The wavefunction in Eq (14) can be obtained from 
Eq. (13) by the unitary transformation U = e~ l(JzQ l 2 . 
So, the physical results should remain unchanged while 
using any of the the two wavefunctions. In the follow¬ 
ing our calculations are always consistent in the limit 
n = 0 (no disclination) with the results obtained in the 

hterature 1 ^®!] 

A. Eigenvalues and Eigenfunctions 

By writing now explicitly H4/ r (p, 6) = FT r (p, 0), we 
obtain 


1 ( <f> n 

s X ( J + ^ + i T 

1 ( <f> n 

St( J+ ^ + I T 


$b = (£-t6)$ a (15) 

= (£ + t5)$b ( 16 ) 


with £ — E/Eq , where Eq = and S = The 


energy spectrum is given by solving Eqs. (15) and (16) 


^n,J — S ]j ^ ^ + q T ) ( 17 ) 


which may be also written as 


J 


77/ l 

£n,J = s\l$ 2 + + p + ! + Jr)(m + /3 + -■r) + ^ 


n 
4 ‘ 


(18) 


where we have defined the reduced magnetic flux /? = 
and the labels s = +1 and s = — 1 which refer to electron 
and hole bands respectively. The energy spectra for K 
and K' valleys are not same. The valley degeneracy (r) 
is broken by the presence of the disclination. This broken 
degeneracy has important effects on charge currents and 
wave packet revival time, which we will discuss later. The 
amplitudes, reminiscent of the radial part of two spinor 


components are 

= 1, (W) 

= FTrf 55 (20) 

where we defined, v T — J+ ^ + 4 r . 

Here we fix the parameters A = 20 meV and p = 
50 nm. Fig [2] shows the energy level (for reasons of clarity, 
in the legends of figures we will denote physical quantities 
with K or K' for the two Dirac points (like in £k) instead 
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FIG. 2. Energy levels (a) 8 k and (b) 8k ' of a single layer 
graphene quantum ring as a function of a mass term 5 with 
P — 0 and m — 0. Blue, green and red curves are for graphene 
with n — 0,1, 2 respectively. 


FIG. 3. Energy levels 8k for a graphene layer with a square 
defect as a function of m for different values of the reduced 
magnetic flux ft. The solid curve is for ft = 0, large dashed 
curve for /5 = 1.0 and small dashed curve for P = —1.0. 


K and K' points respectively. Graphene with a negative 
curvature (n is negative) exhibits the opposite effect. As 
for example, graphene with n = — 1 becomes gapless for 
and = — |ff> 0 at K and K' points respec¬ 
tively. 

In Fig. [3] we present the results for the energy levels as 
a function of the angular momentum m for P = 0, 1.0 and 
— 1.0 respectively. From Eq. (18), the energy spectrum 
exhibits a minimum for m = —+ p + jt) for a given 
value of P and it is independent of S. So, the gauge field 
due to curvature shifts the energy minimum and has an 
opposite effect at K and K' points. 

In Figs © and © we plot the energy as a function of 
the reduced magnetic flux (P) for different values of m. 


of the valley index r) of a graphene ring for one valley, 
r = +1, as a function of the rescaled substrate interaction 
S in absence of magnetic field. From Eq. (18), the energy 


difference between the valence and conduction bands is 
estimated to A 8^ j = ^-(J+P+jt) for the value of S = 
0. So, the energy difference between two levels depends 
on the value of the disclination index n for a given value 
of quantum number J. It also takes different values at K 
and K' points. Interesting features can thus appear. For 
example, at the K point, the graphene sheet with square 
defects is gapped (insulator) but at the other Dirac point 
K' it is gapless (metallic). 

In presence of a finite value of magnetic flux, the sys¬ 
tem becomes gapless for a fixed value of P = — (J -h f r) 
i.e. <f> = — (m + | f r)ff>o (for (5 = 0), otherwise the 
system remains gapped. For example, for m = 0 the 
graphene sheet with a pentagonal defect becomes gap¬ 
less for = — |ff> 0 and = — ^<£>0 at K and K' points 
respectively. Similarly, a graphene layer with a defect of 
index n = 2 becomes gapless for = — ff>o and = 0 at 


B. Pseudo-spin polarization 

The pseudo-spin polariz ation is a contribution to the 
orbital angular moment un J 34|35 i It can be calculated 
from the eigenvectors. The components of the pseudo¬ 
spin polarization are given by 


(<7x) = 0 

K) = 


2£r 


(°z) = 


l + £r 

l~£r 2 

1 + £ t 2 


( 21 ) 


and are such that ( a x ) 2 + (a y ) 2 + (a z ) 2 = 1. In Fig. [6] 
we show the curvature effect on pseudo-spin polarization 
(Gy) of a disclinated graphene. It is seen that at p = 0, 
positive curvature has a larger spin-polarization. The 
blue curve in Fig. [6] is for a monolayer graphene ring. Due 
to the ring geometry it also has a permanent pseudo-spin 
polarization at P = 0. 
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FIG. 4. Energy levels (a) 8 k and (b) Ek' for a graphene layer 
with a square defect as a function of (3 = . The blue curve 

is for m — 0, red for m — —1 and mn — 1, green for m — — 2 
and mn — 2, magenta for mn — — 3 and m — 3, respectively. 


C. Charge Currents 



P 



- 2-10 1 2 
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FIG. 5. Energy levels (a) 8k and (b) 8k> for a pentagonal 
defect in a graphene layer as a function of ft = ^. The color 
significance is the same as in Fig. @ 


The angular velocity is defined according to vq = 
0), 0] = a y . So, the current is given by 

Jo = + Qa^b*) = ( 22 ) 

The expressions of the angular current in valleys K 
and K' are found to be 


and 


t k _ 2 (J + P + f) 

8 ~ w n s K + s 


(23) 


rK > = 2 (J + P-5) 
9 ni s K > - 6 


(24) 


We observe that there are three contributions, first a term 
due to the orbital angular momentum (J), second a term 
due to the presence of the magnetic flux (/3) and finally 
the defect contribution (n). It is worth mentioning that 



p 


FIG. 6. Pseudo-spin polarization ( a y ) as a function of f3 for 
m — 0 at point K. The red, green, blue and magenta curve 
are for graphene with n — 2,1, 0, — 1 respectively. 
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FIG. 7. Charge current Jo for a pentagon graphene ring. 
Fig(a) shows the charge current for the valley r = +1 and 
fig(b) for valley r — — 1. Fig(c) is for the total charge currents. 


FIG. 8. Charge current Jo for a graphene ring with square 
defect for r = +1 and r — — 1 in fig(a) and hg(b) respectively. 
The total charge current shows in fig(c). 


the defect contributions are opposite in the two valleys. 
The total current is now 

j e = jf + jf 

_ 2 (J + P)(£k + £k>) - j[{£k - £k>) - 2<51 

" C 2 (Sk + 6)(£k'-6) 1 j 

In Figs. [7] and [8] we have shown the charge current of a 
graphene layer with defects of indices n = 1 and n = 2. 
Note that the current contribution from K and K r are 
not same for pentagon graphene. This is also true for 
odd n. For even n, the contribution from the K —valley 
and K f are the same and they oscillate in phase, this is 
also true for n = cP. The magnetic flux periodicity of 
the charge current depends on the index n. 


III. REVIVAL AND CLASSICAL PERIODICITY 

The revival phenomenon has attracted broad inter¬ 
est over the past decade. Revivals, and associated frac¬ 
tional revivals, are investigated theoretically and also ob¬ 
served experimentally in Rydberg wave packets in atoms, 
molecules, molecular vibrate states and Bose-Einstein 
condensate^l^. The revival method is used to con¬ 
trol the wave packetP§JH21 ? or for isotope separation 43 45 
to give examples of applications. Revivals happen when 
the wave packet comes back to its initial shape during 
its temporal evolution. The time at which the revival 
occurs is called revival time (Tr). The periodicity of the 
wave packet revivals depends on the energy eigenvalue 
spectrum and is independent of the initial shape of the 
wave function. For an initial wave packet, written as a 
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superposition of eigenstates sharply peaked around some 
level rig, the Taylor expansion of the energy spectrum 
E n > around the the energy E no , is given by 


E" 

E n > ^ E K + E' n ,y - n' 0 ) + -^(n' - ntf + .. .(26) 


where every coefficient of the Taylor expansion in 
Eq. (26) gives an important characteristic time scale for 
the propagating wave packet^ 47 . The temporal evolu¬ 
tion of the localized bound state T for a time independent 
Hamiltonian can be written in terms of the eigenfunctions 
u n > and eigenvalues E n > as 


(27) 

n '=0 


with a n > — (u n f |T). If the coefficient a n > are considered 
to be tightly spread around a large n q \n f — n' 0 \ and 
n' 0 1, and taking into account Eq. (26), one has 


iE n ft/h _ e - i /K E n' Q t + E ' n ' Q i n '-n'o)+E'^/2{n'-n' 0 )‘ 2 + ... ) 

_ e -iuJot-2ni(n'-n'^t/Tci-Zniin'-n' 0 ) 2 t/T R +... 

(28) 


which defines different time scale, that is, Tr = 
47iA/|T T ", | and Tci = 27r h/E f n ,. The propagating wave 
initially evolves quasiclassically and oscillates with a pe¬ 
riod Tec then it will spread and collapse (delocalize). At 
later times, given by integral multiples of Tr/2, the wave 
packet will regain its initial shape and oscillate again with 
a period Tci . Moreover, at times that are rational frac¬ 
tions of Tr, the wave packet splits into a collection of 
scaled and reshifted copies called fractional revivals. A 
longer time scale can be defined beyond Tr, the so called 
super-revival time, at which a new cycle of full and frac¬ 
tional revivals emerges again 4 ^. 

We shall construct the initial wavepacket as a linear 
combination 


*(/£>,#) = yy m $ m (p, 6 o (29) 

rri 

centered around a given eigenvalue T mo , with Gaussian 

/ ^ _ (m-mp) 2 

distributed coefficients c m = w 2 a 2 , where a 

is the variance of the distribution. We can write the 
temporal evolution of the initial wavepacket as 


*t{p,0)=E*m{p,0)t 


, iEmt/Ti 


(30) 


The study of the time dependence of the wavepacket is 
calculated in terms of the autocorrelation function, that 
is, the overlap of the initial state To in Eq. (29) and its 
temporal evolution in Eq. (30): 

A(t) = (T 0 |T t > = / T 0 (p,0)T t (p,0)d0. (31) 

Jo 


If we consider the expansion of Eq. ( [26] ) , A(t) can be 
written as 


A (t) = E\ c m\ 2 e iEmt/n . 


(32) 


The occurrence of revivals corresponds to a complete 
overlap with the initial state, for which \A{t)\ 2 returns 
to its initial value of unity. In this case the classical pe¬ 
riod and the revival time can be calculated from Eq. (26) 
and Eq. (18) 


Tci 


27T pfln r <$ 2 n% 

v F V (m+5 + fr) 2 


and 


T r = 


47T/>f2 2 [<5 2 + UT(m + \ + f r) 2 ]: 


v F 


S 2 


(33) 


(34) 


for ft = 0. The characteristic times are depending on 
the type of disclination and are different at K and K' 
points. In Fig. [9] we show plots of the time dependence 
of \A(t)\ for pentagonal and square defects in graphene. 
In all case we choose p = 50 nm, A = 50 meV and 
mo = 5. For pentagonal graphene the initial wavepacket 
has a fast quasiclassical oscillatory behavior with a pe¬ 
riod Tci = 0.2668 ps; this quasiclassical amplitude de¬ 
creases and later there is regeneration of the quasiclassi¬ 
cal behaviour at half the revival time Tr/2 = 32.78 ps. 
Similarly, for a square defect Tci = 0.2133 ps and 
T r = 89.35 ps (see table |T|. 


TABLE I. Classical and Revival times at the two isotropic 
valleys K and K' for a graphene quantum ring with index 
n — 1 and n — 2 with p — 50 nm and A = 20 meV. The 
initial wavepacket were built as a superposition of Guassian- 
distribouted state of width a — 1.5. 


m 

Tci(K) 
n = 1 

Tci(K'), 
n = 1 

Tci(K) 
n = 2 

Tci(K') 
n = 2 

5 

0.2668 

0.2681 

0.2133 

0.2146 

10 

0.2624 

0.2626 

0.2112 

0.2114 

15 

0.2614 

0.2615 

0.2108 

0.2108 

20 

0.2611 

0.2611 

0.2106 

0.2106 

m 

Tr{K) 
n = 1 

Tr(K'), 
n — 1 

Tr{K) 
n = 2 

Tr(k') 
n = 2 

5 

65.56 

50.68 

89.35 

53.20 

10 

407 

354 

539 

406.71 

15 

1268 

1152 

1650 

1360.94 

20 

2888 

2685 

3721 

3216.25 


It is also seen from the results listed in table H that 
the classical and revival times are different at K and K '. 
This difference in time scales actually destroys the oscil¬ 
lating behaviour of the total charge current as shown in 
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FIG. 9. Time dependence of \A(t) \ for initial wavepacket with 
m — 5 and a = 1.5 at the K point. Upper pannel: Tr — 
65.56 ps with pentagonal defect. Lower pannel: Tr — 89.35 ps 
with a square defect. The difference lies in the time scales 
(horizontal axis). 



FIG. 10. Time evolution of J^ for a pentagonal defect (upper 
pannel) and a square defect (lower pannel). The difference lies 
in the time scales (horizontal axis). 


Fig.[lT]for a pentagonal defect in graphene. This unusual 
behaviour is one of the main results of this paper. 

To investigate the time evolution of angular current Jq 
in Eq. (25), we expand in the base of 0 m (p, 6). One finds 
for the expected temporal evolution of the current!^ 


jf = \ c m,m-i | 2 Sin[(£ m - E m -i)t/h\ (35) 


m= 1 


One can find E m - E m _i ~ E' mo + £^ 0 (ra - ra 0 ) from 
which the scales Tci and Tr are arise. It is shown in Fig. 
(10) the time behaviour of the current Jq for defects of 


indices n = 1 and n = 2. When after a few periods the 
initial wave packet enters the collapse phase and the qua- 
siclassical oscillatory behaviour of the currents vanishes, 
only to emerge later at half the revival time. 


IV. SUMMARY AND CONCLUSIONS 

In this paper, we have shown how a specific type of 
topological defect, a disclination which produces an ef¬ 
fective gauge field coupled to the graphene electrons can 
modify the physical properties of the material by the 
breaking of valley degeneracies. 



0 10 20 30 40 50 60 70 


t(ps) 

FIG. 11. Time evolution of total angular current Jq for a 
pentagon defect. 


• The presence of the disclination (the effect of which 
is measured by a non-vanishing integer index n) 
modifies the energy spectrum (which in particular 
can be gapless at one Dirac point and gapped at the 
other) and has the effect of shifting the magnetic 
flux in opposite directions at the two Dirac points 
K and K'. 

• As a consequence, there is a modification of the 
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charge current, and the fictitious magnetic field 
produced by the disclination having opposite signs 
at the two Dirac points, the currents there par¬ 
tially cancel each other. The total current as a half 
periodicity. Although for even n this scenario is 
different, contribution from two Dirac points are in 
same phase and enhance the total charge currents. 

• The disclination may also enhance or smoothen 
(depending on the sign of n) the pseudo-spin po¬ 
larization which is a contribution to the orbital an¬ 
gular momentum. 

• This also leads to remarkable features on quantum 
oscillations and revival times of wave packets, for 
example the fact that the revival times are not ex¬ 
actly identical at the two Dirac points. As a con¬ 
sequence, the oscillating behavior of the current is 
partially destroyed with time. This is true for any 


type disclinations in graphene. 

There should exist experimental signatures of these 
features. For example in the presence of a single square 
disclination, the graphene sheet is gapped at K —point 
but gapless at K' and thus the conductivity of the metal¬ 
lic state should be half that of the pure graphene sheet. 
The change of periodicity in the current dependence on 
the magnetic flux is even a more obvious signature. 
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